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I. INTRODUCTION 



Recently, the CLEO collaboration ^1,] has measured the 7*7 form factor F^*^'"{Q'^ 

with high precision. This data has been processed by Schmedding and Yakovlev (SY) 
using light-cone QCD sum rules (LCSR), taking also into account the perturbative QCD 
contributions in the next-to-leading order (NLO) approximation. In this way SY obtained 
useful constraints on the shape of the pion distribution amplitude (DA) in terms of confidence 
regions for the Gegenbauer coefficients 02 and 04, the latter being the projection coefficients of 
the pion DA on the corresponding eigenfunctions. Note that SY have extended to the NLO the 
LCSR approach suggested before by Khodjamirian pl for the leading order (LO) light-cone 
sum rule method. 

The present analysis gives further support to the claim, expressed by the above mentioned 
authors, that LCSRs provide the most appropriate basis in describing the form factor of the 
7*7 — > 7r° transition. This is intimately connected with peculiarities of real-photon processes 
in QCD 0, 0|- But the method of the CLEO data processing, adopted in 0, seems to be not 
quite complete from our point of view. We think that an optimal analysis should take into 
account the correct ERBL evolution of the pion DA to the scale Qexf^^^ process (the latter 
not to be fixed at some average point, /zsy = 2.4 GeV, as done in ^) and to re-estimate the 
contribution 5'^ from the next twist term. The influence of both these effects appears to be 
important and it is examined here in detail. Furthermore, we are not satisfled with the error 
estimation performed in the SY analysis, for reasons to be explained later, and prefer therefore 
to use a more traditional treatment to determine the sensitivity to the input parameter 5^ and 
the construction of the 1-cr and 2-a error contours. 

Our main goal in the present work will be to obtain new constraints on the (02, 04) DA 
parameters from the CLEO data, taking into account all the remarks mentioned above, and 
then to compare them with the constraints following from QCD SRs with nonlocal condensates 
(NLC). We will not repeat here the derivation of the main results of LCSRs, as well as those 
related to NLC SRs, but we will refer the interested reader to 0, 0] and correspondingly to 
0, 0] and references therein. But for the sake of convenience we included a technical exposition 
of our approach in comprehensive appendices. 

The paper is organized as follows. In Sec. |n] we review different QCD approaches to 
calculate the transition form factor F"'' "''^ [Q"^) , having recourse to QCD "factorization theo- 
rems" 0], encompassing both perturbation theory and LCSRs. The analysis of the CLEO 
data is discussed in Sec. Illll in conjunction with the SY approach in comparison with other 
approaches/approximations. In Sec. IIVI we present a complete NLO analysis of the CLEO 
data with a short discussion of the ELM setting procedure. Sec. El includes a comparison 
of the QCD SR pion DA models with the results obtained in Sec. IIVI from the CLEO data 
processing. In Sec. IVII we summarize our conclusions. The paper ends with flve appendices: 
in Appendix A we re-estimate the value of the twist-four scale 5^. In Appendix B the old 
Chernyak-Zhitnitsky (CZ) result for 02 is discussed, paying attention to evolution effects. In 
Appendices C and D the two-loop results [1, |^ for the purely perturbative part of the form- 
factor calculations and the ERBL evolution of the pion DA are outlined. Finally, in Appendix 
E, all needed calculation details for the NLO LCSR are presented. 
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II. TRANSITION FORM FACTOR F^*^'^(Q2) WITH LCSR 



A. Factorization of the F''^ ''^ form factor. Standard results 

The form factor of the process 7*(q'i)7*(5'2) 7r°(p) is defined by the matrix element 



(2.1) 



where gi + g2 = = —qi > with = — > are the virtuahties of the photons, 

i/i = i^u-y^u — ^d'y^d) is the quark electromagnetic current, and vr°(p) is the pion state with 
momentum p. If both virtuahties, and g^, are sufficiently large, the T-product of the 
currents in (j2.ip can be expanded near the light-cone z"^ = 0. This expansion results in 
factorization theorems, which control the structure of the form factor at leading twist. As 
a result, the form factor can be cast in the form of a convolution over the momentum fraction 
variable x (of the total momentum p) to read 



^QCD''(<5^ = U jdx T{Q'\ g^; /ip; x)ip^{x] /ip) = fnT{Q'^, g^ /ip; x) <S) '^nix; /ip) . (2.2) 



The hard amplitude T[Q^, g^; x) - playing the role of the Wilson coefficient in the OPE - is 
calculable within QCD perturbation theory (PT): 



T = To + -J^T, + 
An 



«s(/iR) 

Air 



T2 + ... 



(2.3) 



while the pion DA (Pt^{x; fip) contains the long-distance effects and demands the application 
of nonperturbative methods. Due to factorization theorems, DAs enter as the central input of 
various QCD calculations of hard exclusive processes. The LO term T^i in ()2.3p depends only 
on kinematical variables, but the NLO amphtude Ti, calculated in j^l^], depends also on the 
factorization scale 



To{Q',q';x) = NtCo{Q\ q' ; x) = Nt 



Q2j. _|_ q2(^^ _ r^^ 



ix ^ 1 — x) 



Ti(Q^g^;/i^;x) = Nt 



In 



2 ; ■Co(g',g';u)®K)(M,a;) + Ci(g",g";x) 



(2.4) 
(2.5) 



Here, Nt = {el - ej)/V2 = V2/Q is the QCD normalization factor, = - [(gi - g2)/2]^ = 
(Q^ + Q'^)/2, and and denote, respectively, the scale of renormalization of the theory 
and the factorization scale of the process. Co,i and Vo,i are, respectively, the perturbative 
expansion elements of the coefficient function C of the process and the ERBL kernel V in the 
LO (subscript 0) and NLO (subscript 1) approximation ^. The structure of Ti is discussed 
in more detail in Appendix Ol- Here we only recall those features relevant 

for our discussion: Eq. ()2.5p . which specifies the definition of the coefficient function Ci 
is a consequence of the QCD factorization theorems 0|. Due to these theorems, the first 
logarithmic term in Eq. (j2.5|) . originating from coUinear divergences, can be absorbed into the 



wherein Vi is needed to evolve the pion DA in the NLO approximation. 
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renormalization of the DA, (pTr{x]fip), following the ERBL equation (see App. iDl Eq. ()D.1|) ). 
Namely, the formal solution of the ERBL equation in the 1-loop approximation is 



/ 2\ ERBL RC/ ■2\ 



exp 



(2.6) 



(here /?(a) = —boo^ is the 1-loop /3-function)^. We adopt here as a factorization scale /i| = Q^. 
For that case, F^*^*'^ , obtained in Eq. ()2.2j) at the scale /i^, can be transformed into 



^7 7 7r 



^RG 



To{Q\q';x) + ^^T,{Q\q';Q';x) 



V.^«(x;Q^). (2.7) 



To fix the renormalization scale fi"^, one needs to go beyond the NLO approximation. In the 
absence of such information and in order to further simplify the NLO analysis, we also set 
f^R — f^F — Q^- is useful to expand ip^^{x; /x^) over the eigenf unctions '?/'„(x) of the one-loop 
ERBL equation, i.e., in terms of the Gegenbauer polynomials C^^^(O) 



n=0,2,4,... 



ao 



1; iPnix) = 6x(l -x)C3/2(2x 



(2.8) 



In this representation, all the dependence on /i^ is contained in the coefficients an{^^) and is 
dictated by the ERBL equation. Different reasons, explained in and [Si, |6||], point to the 
possibility of retaining only the first 3 terms in this expansion. Following then this approx- 
imation, the form factor can be parameterized by only two variables ci2(/^o) ^^"^ '^4(/^o) ^^^^ 
accumulate the mesonic large-distance effects (at some scale /Xq)- 

A special case of the OPE appears when one of the photons (^2) is near the mass shell. 
It is instructive to present here an evident NLO perturbative expression for Eq. ()2.7|) in the 
formal limit ^ (cf. Eq. ()C.8|) in Appendix [n)). 



F^^'^'m = {1 + S(Q^) + ^4f^Cp [-5 + 7.9 • S(g2) - 3.8 ■ A(Q^ 



80 5^ 

■27~g2- 



(2.9) 



with S(Q2) = a2(Q2) + a4(Q2) and A(g2) = a^iQ^) - a^iQ^). 

On the r.h.s. of Eq. (|2.9|) . the twist-four contribution is included in terms of the twist-four 
asymptotic DAs, presented in [3j]. The coefficients 02,4 (Q^) and the twist-four coupling ^^(Q^) 
are evolved to the scale using, correspondingly, the renormalization group (RG) equation 
to NLO and LO accuracy^. 



B. Why do we need LCSRs for the transition form factor? 

A straightforward calculation of F^*"''^{Q'^) in QCD is not possible. In particular, at the 
formal limit ^ 0, it is not sufficient to retain only a few terms of the light-cone OPE 



^ and one should mean 



(g) Vo('u„_i, u„) (g) 9?Tr(u„; 



The reasons for this treatment will be discussed in more detail in Sec. |l3 
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of ()2.H) . One has, in addition, to take into account the interaction of the small- virtuahty 
photon at large distances proportional to 0(1/ \/q^), (for a recent discussion, consider section 
4 in |3] and [l^). The LCSR method allows one to avoid the problem of the photon long- 
distance interaction by providing the means of performing all QCD calculations at sufficiently 
large and then use a dispersion relation to return to the mass-shell photon. To isolate the 
dangerous neighborhood of = 0, one should apply an appropriate realistic model for the 
spectral density at low s, based, for instance, on quark-hadron duality. Using this sort of 
analysis and by employing analyticity and duality arguments, the following expression was 
obtained in 



T 



LCSR 



ds ^ 
— Im 

nip 



„ -^QCD ' ^) 



_ -^QCD W 5 



(2.10) 



so 



F^nr-^ {Q"^ 1 s) on the r.h.s. of ()2.10|) is taken from Eq. ()2.2|1 and the Borel parameter is 



QCD 



0.7 GeV^, whereas 



1.5 GeV^ denotes the effective 



p is the p-meson mass and sq 
threshold in the p-meson channel. 

This program has first been suggested in general form in and was realized there for the 
LO approximation of the process. Taking -Fqc^'^ at LO in a^, the form factor -Flolc('5^ 
obtained from Eq. (j2.10|) with 



was 



77.7 7 T 
-^LO QCD 



fnN: 



T 



Co{Q'^,q'^;u)(g)(p^{u) 



[CoiQ 



2 2 
q \u 



and 



Im 



TT 



77-7 7 -n- 
-^LO QCD 



fnV2 



S + Q2 Q2 ds 



(2.11) 



(2.12) 



u=Q2/{s+Q2) 



Here, the numerically important twist-four contribution was also included using a simple 
asymptotic expression for the twist-four DA contribution 



(2.13) 



Note that the coupling 6"^ is determined from the matrix element 

{7rip)\gjG^^ru\0) = iS^Up^ . (2.14) 
The estimates for ^^(yU^), including also its RG-evolution, are presented in Appendix 1X1 



C. The framework of NLO LCSR for the form factor 

An application of the above scheme to the NLO was more recently performed in The 
tt.^-corrections to this process are expected to be rather large, of the order of 20% (see, e.g., 
[10|). The size of these corrections can be roughly estimated from the perturbative expressions 

Q/ ____ 

for their asymptotic parts by comparing to each other the values 1 and — 5— Cp in Eq. (12.91) . 

Therefore, for a quantitative description of the form factor, this contribution should be taken 
into account. This important step was done by Schmedding and Yakovlev, who have used the 
NLO perturbative expression for -Fqcd'' ^'i- (EUDl), 

F^^olcBiQ^ /^') = F^o'ocBiQ^ q') + ^T^^Tl(Q^ q'; /i^; x) ® ^.(x; p^) , (2.15) 
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where fi"^ = fip = fi^ ^ Q^, to construct a NLO version of the LCSR for the form factor 
pi*i'^(^Q'^y The spectral density Im[Fj^Lo qcd('5^' ■^j Z^^)]/'""' based on Eq. (|2.15p . depends on 
the scale /x^. In the original paper 2], this scale was fixed by relating it to the mean value of 
with respect to the CLEO experimental data, i.e., by setting jj? = /igy = (Qexp) ~ (2.4)^ GeV^. 
Use of Eq. fl2.15|) implies that one accounts for the scale dependence in ip-jrix] Q^) (e.g., for the 
different CLEO points) via the leading order perturbative formula 

rather, than using the RG expression, given by Eq. ()2.6|) . This seems to be a rather crude 
approximation, given that other reference points are evolved with fi^ by utilizing the NLO 
ERBL evolution equation. 

III. CLEO DATA ANALYSIS REVISITED 

Let us scrutinize the form factor approaches, discussed above, in close comparison with the 
CLEO data Q. 



1 + 



47r 



■In (QV/i2) .\/o(x,m) 



A. Theoretical approaches to the CLEO data 

In this subsection - in order to be in close analogy with the original SY paper - we shall 
adopt their scale definition of the CZ DA. It is worth noting, however, that the genuine CZ 
DA differs from that definition because it is determined at the much lower normalization scale 
yUo = 0.5 GeV (a discussion of this important point is relegated to AppendixlHl). To distinguish 
the two models in the present analysis, we will use in what follows a special notation: (i) CZ 
DA parameters originating from the SY prescription ^(/i^ = 0.5 GeV^) = 2/3), and by using 
a 2-loop evolution to the scale /igy, will be denoted with the superscripts CZSY, whereas (ii) 
those conforming with our prescription ^(/U^ = 1 GeV^) = 0.56, see Appendix^ and being 
2-loop evolved to the scale /Xgy will be marked by the superscripts CZ. 

The status of these approximations (NLO PT, LO LCSR and NLO LCSR) is illustrated in 
Fig. n] by the relative positions of the corresponding admissible regions for these parameters in 
the (02, 04) plane. Here, the regions enclosed by the needle-like- and ellipse-like solid contours 



correspond to a la-deviation criterion (CL=68%) while the broken contours refer to 
a 2(T-deviation criterion (CL=95%). Note that these contours have been derived by taking 
into account only the statistical error bars in ^ (see their Table 1). This marks a crucial 
difference between our processing of the data and that in j3| , where the "theoretical-systematic 
uncertainties" have been involved in the statistical analysis together with the statistical ones. 
In other words, we do not "smear" the quantities 6^ (or = 1) over their corresponding 
(theoretical) error-bar intervals. In our opinion, such a manner would require an additional 
substantiation and further suggestions about the distribution of these errors that we want to 
avoid in our analysis. Hence, instead of that, we process the data at a few fixed values of 5^ 
in order to clarify the sensitivity of the results to this parameter. 

It should be clear that a really admissible region might be somewhat larger than the pre- 
sented "purely statistical" contours in Fig. [TJ a price one has to pay for our strict way of data 
processing. 
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FIG. 1: The plane (02, 04) and results of different analyses of the CLEO data for the form factor F^^^(Q^). 
We show here the best-fit points of LO LCSR, NLO LC SR (SY) and a NLO PT analysis for two values of the 
twist-four coupling = 0.2 GeV^ (a) and = 0.18 GeV^ (b). The solid contours correspond to a Icr-criterion, 
while the broken ones to a 2CT-criterion. We also plot for comparison the positions of the CZSY (■) and the 
asymptotic (♦) DAs. All values are evaluated at = (2.4 GeV)^. The inclined straight line on the right 
hand-side of the figure is the modified diagonal 02 -|- 04/8 =C0NST, corresponding to the new la-ellipse of the 
NLO LCSR analysis. 



(i) The large contour, enclosing the origin, corresponds to the LO LCSR form factor, dis- 
cussed in 0. Note that the point denoting the position of the asymptotic DA lies just 
on the periphery of the Icr-contour and that the CZ DA one (see [3j), lies far outside. 

(ii) The needle-like contour on the left top corner of both parts of Fig. ^ corresponds to 
Eq. ()2.9p and is stretched along the "diagonal" 02 + 0.75-04 = CONST. The weak 
dependence of Eq. fl2.9|) on A = 02 — 04 slightly turns the angle of the diagonal ~ 
371/4 — > O.Svr. On the other hand, taking into account the evolution with of 02 
and 04 for every Qg^p inakes the contour finite - much like a "diagonal" needle-like 
strip. As we have mentioned above, the formal limit of the NLO expression, Eq. (j2.9p . 
cannot give a reliable result for the form factor. In this context it is interesting to 
mention that the corresponding contour is located outside the regions determined by the 
LCSRs. Furthermore, all known DA models (see, e.g., Table|l} and the phenomenological 
predictions (0,01 and references therein) are located far away from this contour - clearly 
demonstrating the poor reliability of the corresponding perturbative approach. 

(iii) At least the heavy-line contours (enclosing also the SY point 0) correspond to the SY 
approximation. These contours do not overlap with those corresponding to the LO LCSR 
ones - even at the 2(T-deviation level. Therefore, ag-corrections are crucially important 
in extracting the DA parameters. Our best-fit point with respect to the NLO LCSR is 
close to but not coinciding with the one presented by SY (compare entries 3 and 4 in 
TableU}. In Fig. ^a) the full circle inside the contour is just the SY point. The best-fit 
points, corresponding to different approximations and models, considered in the present 
analysis, are collected in Table HJ where the notation xfdoi) ~ has been used (in 
correspondence to the number 15 of the CLEO experimental data points). 

It should be stressed that the admissible region (heavy-line contours), obtained with our data- 
processing procedure, differs from that in Our contours look slightly larger than theirs 
despite the fact that possible theoretical/systematic uncertainties were not included in our 
consideration. Just because of this latter reason, our contours possess another orientation 
relative to those of SY, as one appreciates by comparing Fig. (TJ^a) and Fig. |21 with Fig. 6 in 0. 
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TABLE I: Compilation of best-fit points shown in Fig.IHa) with 5'^{1 GeV^) = 0.20 GeV^. 



Best-fit points/models 


02 (^Sy) 


«4(^iY) 


v2 

A(d.0.f.) 


NLO PT (irni) best fit 


-0.44 


+0.40 


0.49 


LO LCSR best fit 


+0.04 


-0.19 


0.48 


NLO LCSR best fit 


+0.20 


-0.17 


0.48 


SY LCSR [2J 


+0.19 


-0.14 


0.49 


BMS model \5] 


+0.13 


-0.08 


0.74 


CZSY DA \2} 


+0.28 


-0.009 


1.8 


CZ DA [1^ 


+0.35 


-0.006 


3.6 


Asympt. DA 


-0.006 


+0.00 


3.0 



Moreover, the CZSY DA model appears to be seemingly closer {x^ = 1.8) to the best-fit point 
(0.20, —0.17) than the asymptotic one (x^ ^ 3). But the genuine CZ DA with a^^i^ily) = 0.35 
and a4^(/igY) = —0.006 (consult the discussion in Appendix |B|) generates a value of = 3.6, 
which is larger than that of the asymptotic DA. 

The best-fitted linear combination^ of a2, a^, that determines the large axis of the NLO 
LCSR contour (see Fig. HJ^a)) and parameterizes its orientation, is found to be 

as + - ■ 04 = 0.143, (3.1) 
3 

instead of 02 + 0.6 ■ 04 = 0.11 ± 0.03, reported in 0. Note that the SY point also belongs to 
the diagonal: af^ + 1/3 ■ af^ ^ 0.143. The coefficient 1/3 in (jSH) can be predicted without 
any fitting; it is solely determined by the structure of the NLO LCSR (|2.10p . Indeed, Eq. 
fl2.10|l can be rewritten as 

AoiQ^f^lr) + MQ^f^lr) ■ «2(/xiY) + A,{Q\fily) ■ a^if^l^) = ■ F^''^^Q') 

(to be compared with the fit given in Eq. (|3.1|) ) and the discussed coefficient expresses the 
average value of the ratio A^^Q'^ , fi'^y) / A2{Q^ , f^sy) ■ Notice that this ratio, averaged over the 
CLEO data range {Ql^p}, amounts to 0.31, while the r.h.s. of Eq. ()3.1|) is determined by the 
experimental data on the form factor. The coefficient 0.6, obtained in the SY fit 0, can be 
associated with the same ratio at the mean value /Xgy = (Qexp)- The ratio A^i^Q'^) / A2{Q'^) is 
a concave function in and, therefore, its mean value, {A^IQ"^) / A2{Q'^)) , is smaller than its 
value, A^^fi'^y) / A2{figY) , at the "mean point" /igy. 

B. Sensitivity to input parameters 

As it turns out, the location of the admissible 02, regions is rather sensitive to the value 
of the input parameter 5^. To illustrate this point, we have repeated the data processing with 
an admissible (near its low boundary) value 5^ = 0.18 GeV^ (see Appendix E}- AH contours in 



^ Dubbed "diagonal" in what follows. 
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Fig.^b) shift closer to the asymptotic point (♦), but their relative positions do not drastically 
change and, hence, the main conclusions (i-iii) of the previous subsection remain valid. The 
results of this data processing are presented in Fig.^b) and in Table ITTl One appreciates that 
the hierarchy of the different models (lower parts of both Tables) with respect to the NLO 
LCSR best fit does not change, though the values of can change significantly. Indeed, the 
point marking the BMS model moves from the 2cT-deviation level at 5^ = 0.2 GeV^ (see Table 
U} inside the la-deviation region near the SY LCSR point at = 0.18 GeV^ (cf. row 5 in 
Table ITT]) . Therefore, the value of 5^ and, in general, also the value of the twist-four term can 
substantially affect the locations of the admissible regions. But all other options, like the CZ 
DA and the asymptotic DA, remain excluded at the 2(T-deviation level. 

TABLE II: Compilation of best fit points shown in Fig. Ctb) with 5^(1 GeV^) = 0.18 GeV^. 



Best-fit points/models 


«2(Msy) 


«4(^Sy) 


v2 

A(d.o.f.) 


NLO PT ((Tll) best fit 


-0.37 


+0.34 


0.49 


LO LCSR best fit 


-F0.006 


-0.14 


0.49 


NLO LCSR best fit 


+0.17 


-0.14 


0.48 


SY LCSR \2] 


+0.19 


-0.14 


0.51 


BMS model \5] 


+0.13 


-0.08 


0.56 


CZSY DA [2j 


0.28 


-0.009 


2.2 


CZ DA [1^] 


-hO.35 


-0.006 


4.2 


Asympt. DA 


-0.006 


-FO.OO 


2.3 



Let us pause for a moment and turn our attention to a recent paper by Diehl et al. jl^ . The 
authors of this work employ a purely perturbative QCD approach to analyze the CLEO data 
without taking into account the twist-four contribution, i.e., using Eq. ()2.9|) with ^^(Q^) = 0. 
They consider this treatment justified given the possible large uncertainties in estimating the 
twist-four contribution (which in the SY procedure is taken to be ±20%). Comparing their 
results with those of Schmedding and Yakovlev [3], Diehl et al. correctly note that the relative 
weights of (32 (/^o) '^4(/^o) F"'*"''^ display a much stronger Q^- dependence than in the 
leading-twist case with the consequence that the allowed SY parameter region becomes much 
smaller than in their approach. However, the size of the twist-four contribution is crucial for 
accurately extracting the parameters 02 and 04 - as we have just demonstrated. Therefore in 
our analysis we use a different approach: the value of 5^ is connected with the parameter of 
the vacuum non-locality. We first fix the value of and then we allow for the parameter 5^ to 
vary in a 10% range. The whole uncertainty in 5"^ for the selected range of A^ = 0.4 — 0.6 GeV^ 
amounts then to about 30% in accordance with . This strategy enables us to use the CLEO 
data as a direct measure (a vacuum detector) to select that model for the QCD vacuum, which 
provides the best agreement between theory and experiment. 

IV. COMPLETE TWO-LOOP ANALYSIS OF THE CLEO DATA 

In the previous section we have demonstrated the high sensitivity of the DA parameters 
(02, 04) to strong radiative corrections for the form factor, as well as to the scale of the twist- 
four contribution (see Fig. ^ a(b) and j^). Therefore, to obtain these parameters from the 
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CLEO data in a reliable way, one should take into account the radiative corrections in the 
most accurate possible way. To this end, we want to improve in this section the accuracy of 
the extraction procedure of (02, 04) at the NLO level. A new estimate for 5^, the magnitude 
of the twist-four contribution, is also introduced in the present analysis (see below). We also 
briefly discuss an attempt to go beyond the level of the NLO, having recourse to a recent 
calculation |0| of the radiative correction based on the BLM scale setting. 

A. Complete NLO analysis 

Here we use the complete 2-loop expression for the form factor F'^ ^{Q'^i given by Eq. 
(12. 7j) . For this reason, we put /i^ = in (j2.15|) so that for the quantities 

the NLO evolution is implied. Then, we substitute the spectral density p{Q'^,s;^'^ = Q"^), 
derived in (see the text below Eq. ^J^), in LCSR (imUl to obtain Ft*t^(Q2) in a regular 
manner and to fit the CLEO data over Q"^ G {Ql^p}- The evolution ipj^^x; Q"^) = U{h\y ~^ 
Q^)(y9^(x; /igy) is performed for every point Qcxpj with the aim to return to the normalization 
scale /igY and to extract the DA parameters (02, 04) at this reference scale. Stated differently, 
for every measurement, {Qexp; ^^'''^''i.Qlxp)}^ its own factorization/renormalization scheme has 
been used so that the NLO radiative corrections are taken into account in a complete way. 

The accuracy of the procedure is, nevertheless, still limited owing to the mixing of the NLO 
and the LO approximations. Indeed, the value of the twist-four coupling 5^(/i^), as well as 
its RG-evolution with /x^, are estimated in the LO approximation. This quantity enters the 
LCSR formula (j2.15p together with the NLO-part and can lead to an additional uncertainty. 
In order to improve the theoretical accuracy of the values of 02 and 04, extracted from the 
CLEO data, one has to re-estimate the twist-four contribution, 5^(1 GeV^), with a better 
accuracy. 

To summarize, our data processing procedure differs from the SY one in the following 
points: 

1. as{^J^) is the exact solution of the 2-loop RG equation, rather than the approximate 
(but popular in the HEP community) expression ()C.15|) [nj, that was used in the SY 
analysis. 

2. All logarithms ln(Q^//i^) are absorbed into the evolution of the pion DA, performed 
separately at each experimental point Qexp (compare Eqs. (j2.6p . (j2.7|) with Eq. (|2.15|) ). 
The corresponding expressions for the form factor are collected in Appendix E. 

3. The value of the parameter 5^ has been re-estimated to read 5^(lGeV^) = 0.19 ± 
0.02 GeV^ (see Appendix A), and this value has been used in the data processing. 

This processing of the CLEO data produces the admissible regions, one of which, corresponding 
to 5"^ = 0.19 GeV^, is shown in Fig. |2fa), where the original SY regions (Fig. 6 in j^) are 
also presented in Fig. I21d) for the ease of comparison. To produce the complete 2a- and 
Icr-contours, corresponding to 5^(lGeV^) = 0.19 ± 0.02 GeV^, we need to unite three regions 
obtained for different values of the twist-four parameter: 5^ = 0.17, 019, 0.21 GeV^. This 
procedure is illustrated in Fig. Efb) using as an example the 2(T-contour. Let us remind the 
reader in this context that the SY contours are stretched along the "LO perturbative" diagonal 
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FIG. 2: Results of the analysis of the CLEO data for the form factor Ft^^*^{Q'^) in terms of fiducial areas 
in the plane (02, a4). (a) We show the best fit point of the new NLO LC SRs (+, termed "New", which is 
associated with the value of the twist-four coupling ~ 0.19 GeV^, and display by the heavy (broken) contour 
the regions corresponding to a Icr (2cr)-deviation. We also show the SY point (•), the BMS point (X), the 
asymptotic point (♦), and the CZ point (■) - all evolved to the SY scale. The inclined straight dotted line 
in the left figure is the modified diagonal 02 + 0.804 = 0.159, whereas the inclined straight dashed line is the 
original diagonal 02 + 04 =C0NST. (b) We show three 2o'-regions (dashed contours), obtained for different 
values of the twist-four parameter: 5^ — 0.17, 019, 0.21 GeV^ and also their unification: 2(T-contour (thick 
line). All points and straight lines are as in Fig.[5fa), except which denotes a brand-new result obtained on 
a transverse lattice 0| . (c) The contours of the previous NLO analysis a la SY and the "old" best-fit point 
(+) are presented for comparison (cf. the lower right corner of Fig. QJa) and use the scale of Fig. [2Ia)). (d) 
The original SY plot (note the different scales on the axes and the symbol ■ denoting the CZSY DA) is shown 
for comparison. All values shown are evaluated at the scale /i — 2.4 GeV. 

a2 + 04 =C0NST (the dashed straight hne on the l.h.s. resembles exactly this diagonal) while 
the solid (dotted) contours correspond to the 1 a (2 a) regions. This stretching of the contours 
appears here because of the SY manner of the data processing, namely, because the theoretical 
uncertainties of the input parameters were also involved in the statistical analysis. The new 
best-fit point (+, "New"), as well as the whole a-contours themselves appear to be displaced 
in Fig. I^a) (approximately) along the new diagonal (cf. Eq. (jSH])), 

a2 + 0.3 -04 = 0.16 ±0.007. (4.1) 

In Fig. Efc) we present for comparison the contours of the previous NLO analysis in the sense 
of SY (low right corner of Fig. ^a)) drawn, however, at the scale of Fig. |2ta). The positions 
of the best-fit points and models are provided in Table IIIIl 
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TABLE III: Compilation of best-fit points shown in Fig. I^a) with d'^{l GeV^) = 0.19 GeV^. 



I— < /~1C i~ Tn i~ T "\ T~l T~ / T^"l 1 o 

OtJoL-ilL pUlllL / lliUQtJlo 






v2 

A{d.o.f.) 


U2 ~r U.o • tl4 


New NLO LCSR best fit 


+0.23 


-0.22 


0.47 


0.159 


SY NLO LCSR \2] 


+0.19 


-0.14 


0.59 


0.153 


BMS model [5] 


+0.14 


-0.09 


1.0 


0.113 


Asymptotic model 


-0.002 


+0.00 


4.3 


-0.00 


CZ model fl3] 


+0.43 


-0.003 


4.3 


0.43 



It should be clear from our discussion that these new contours are somewhat smaller than 
the previous ones (Fig. I2tc)), but slightly larger than the original SY ones (Fig. Efd)), and 
show another orientation along the diagonal Eq. 1)4.11) . The difference between the new regions, 
determined in the present analysis, and those of the SY one is remarkable. For instance, the 
SY point appears now near the boundary and inside the Icr-region in Fig.|21^a). Moreover, the 
preliminary (i.e., for 04 = 0) SY best-fit point, 05 = 0.12 + 0.03 and the phenomenological 
estimates for (02, a^), presented in |ia], (04 = 0, 02 (1 GeV^) = 0.1 ±0.1), lie on the boundary 
of the united 2(T- region (see Fig. I2tb)). 

B. Beyond the NLO approximation: effects from BLM scale-setting 

The renormalization scale fi^ in Eq. ()2.7|) can be fixed by a NNLO calculation of T2 following 
the BLM prescription. Recently, the NNLO contribution to C2, proportional to 60 and required 
for the BLM scale setting, was obtained in for a kinematics with q2 = 0. As an exercise, 
let us perform the new fit to obtain the scale /x^ = /^Ilm and the best-fit point {af a^^^} 
for the NLO expression given by Eq. ()2.9p . We follow the same procedure as in Sec. 1111 Al 
replacing this time ^^(Q^) — > o:s{fi'^i^^{a2, 04, Q^)) in Eq. ()2.9p . where for /Xblm('^25 0,4, Q"^) Eq- 
(7.7) from ^| is used. As it turns out, practically for all points in the considered domain 
in the lower half-plane 04 < 0, the BLM setting leads to the condition /xIlm ^ Q^, in 
conformance with the results of [16| . Therefore, for this region, the BLM setting seems to rule 
out predictions from the NLO perturbation theory. 

Only for points within a rather thin strip in the upper half-plane 04 > (cf. Eqs. (7.2a), 
(7.7) in j3l), the BLM setting gives /u|lm > Q^- Interestingly, the case discussed in Sec. IIII Al 
(ii) and based on Eq. ()2.9|) belongs just to this thin strip. The corresponding (02,04) values 
are displayed below for comparison together with the initial result (second row) without the 
BLM setting. 

«f ""(/iiY) = -0.55, a?LM(^2^) = 0.48 ^^l = /x^m ^ QVO.35 (4.2) 
a2(AiiY) = -0.44, a4(/iiY) = 0.40 fil = Q^. 

To calculate the imaginary part IiyiFq^y)^ {Q"^ , s) , used in Eq. ()2.10p . one needs to know the 
NNLO contribution proportional to bo for g| 7^ 0, which is still not computed. For this reason, 
the results obtained with the BLM scale setting fall out of the region of the NLO LCSR 
analysis. The calculation of the complete NNLO contribution or, at least, its convolution with 
(fias is a very demanding task that has not been accomplished yet. 
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V. PION DA FROM QCD SR VS CLEO DATA 



Let us now turn to the important topic of whether the CLEO data is consistent with 
the non-local QCD SR results for (/j^. We present in Fig. El the results of the data analysis 
for three central values of the coupling S\fi'^ = 1 GeV^) = 0.19,0.235,0.29 GeV^ which 
in turn correspond to three admissible values of the vacuum non-locality parameter = 
0.4,0.5,0.6 GeV^. For each value of from this ensemble, we define the corresponding 
central value of 5^ and its uncertainty (for details, see Appendix^. Then, we process the 
CLEO data as described in the previous section and obtain the complete 2a- and Icr-contours 
on the plane (02,04), following from the CLEO experiment. An example of these regions is 
represented in Fig. I^b) and is also displayed in Fig. Ol^a), where the 2cr-contour is shown as 
a solid line and the Icr-contour dashed one. 

The task now is to compare these new constraints with those following from the QCD SRs 
with nonlocal condensates. We have established in ||5j] that a two-parameter model fnix] (12, 04) 
really enables one to fit all the moment constraints for {C,'^)n that result from NLC QCD SRs 
(see for more details). The only parameter entering the NLC SRs is the correlation scale 
Ag in the QCD vacuum, known from nonperturbative calculations and lattice simulations (for 
a discussion and references, see Appendix . 




0.1 0.2 0.3 0.4 




0.1 0.2 0.3 0.4 




0.1 0.2 0.3 0.4 



FIG. 3: Three 2(T-contours of the admissible regions following from the analysis of the CLEO data for 

different values of 5"^: (a) - for = 0.4 GeV^ and 6"^ = (0.19 ± 0.02) GeV^; (b) - for = 0.5 GeV^ 

and ^2 = (0.235 ± 0.025) GeV^; (c) - for A^ = 0.6 GeV^ and = (0.29 ± 0.03) GeV^ Solid lines in all 

figures enclose the 2(T-contours, whereas the lu-contours are enclosed by dashed lines. The three slanted and 

n 

shaded rectangles represent the constraints on (02, a4) posed by the QCD SRs 5] for corresponding values of 
A^ = 0.4, 0.5, 0.6 GeV^ (from left to right). All values are evaluated at fi — 2.4 GeV after evolution. The 
marked points are explained in Fig. |21 except for the point marked with X. This point represents here the 
BMS solution, which corresponds to a particular value of A^. 

The three slanted and shaded rectangles in Fig. |21 are the constraints on the Gegen- 
bauer coefficients (02, a^j resulting from the NLC QCD SRs at different values of A^ = 
0.4, 0.5, 0.6 GeV^ , 0, |^. The overlap of the displayed regions in Fig. El can serve as a 
means of determining the appropriate value of A^. In fact, one may conclude that the value 
Ag = 0.4 GeV^ is more preferable relative to the higher values of A^. It should be noted, how- 
ever, that even for this lowest value of that scale, the agreement with the constraints in Fig. 
01 is rather moderate and of similar quality as using the SY constraints It is tempting 

to test even smaller values of A^ than 0.4 GeV^ in the NLC SR as an attempt to improve the 
agreement with the CLEO constraints in Fig. El But such values appear to be at the lower 
limits for the A^ estimates from non-perturbative approaches (see Appendix A). Furthermore, 
the NLC SR becomes unstable at such low values of A^. As a result, the accuracy of the DA 
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moments is rather poor and the final constraint on (02, 04) becomes unrehable. Taking into 
account all these arguments, we think that an improvement of the ingredients of the NLC 
ansatz may provide a better agreement with the CLEO data than just using < 0.4 GeV^ . 

VI. CONCLUSIONS 

In this paper we have studied the theoretical predictions for the pion transition form factor 
^7*77r^Q2^ in comparison with the CLEO experimental data ^ on this form factor. We have 
presented a full analysis of this data and contrasted the results with those found in the context 
of QCD LCSR at the NLO level. In this way, we have revised and improved the procedure 
of analyzing the CLEO data, first performed by Schmedding and Yakovlev in [2]. The main 
goal has been to obtain constraints on the shape of the pion DA of twist-2, ipn{x; 02, a^) in the 
most accurate way. The values of the crucial parameters, viz. the twist-four coupling 
and involved in this procedure, have been treated more accurately than in previous 

approaches. The main findings may be summarized as follows. 

1. We have tested different kinds of approximations to calculate F"'*"''^ [Q'^) and revealed 
how the size of and the twist-four corrections can affect the admissible regions of the 
DA parameters (02,04). 

2. New admissible regions for the (02? ^^4) parameters, see Figs. |21I^b), |2^a), - different from 
those in (2] ^ have been obtained, though the constraints do not change drastically in the 
sense that the initial SY best-fit point still belongs to a la-deviation region (CL=68%) 
in this space, whereas the CZ DA and also the asymptotic one are definitely excluded 
at the level of a 2(T-deviation criterion (CL=95%). Moreover, one may appreciate by 
comparing Figs. El^a, b) with Fig. |21^d) that this exclusion with respect to the asymptotic 
DA becomes even more pronounced using our data processing. 

3. The bunch of admissible pion DAs, (/9^*^'^(a;; 02, 04), corresponding to the estimate 
Aq = 0.4 GeV^, that was constructed within the framework of QCD SRs with non- 
local condensates in j^, compares well (at the 2(T-level) with the new more restrictive 
constraints obtained in the present investigation as Fig. El demonstrates. In addition, 
half of the calculated admissible region intersects with the la domain as well, see Fig. 
Ha). 
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APPENDIX A: REVISION OF THE QCD SR RESULTS FOR 5'^ 



The coupling was originally estimated in fill and found to be ^^(^^s _ ^ Qgy2^) _ 

0.2 ± 0.02 GeV^ Here, we re-analyze the QCD SR for 5^, derived in [i^, which is based on a 
non-diagonal correlator of the quark-gluon and quark (pseudoscalar) currents. This SR relates 
5^ to and determines the value of the ratio r = Ag/(25^). Evaluating the SR leads to the 
estimate r > 1 and consequently to Ag/2 > 5^. Moreover, r is rather sensitive to the size of 
the radiative corrections. In this work, we use Ag^ = 312 MeV, obtained recently in a DIS fit 
of the CCFR data in |21i] that leads to ^^^(l GeV^) 0.59. The same sort of analysis in the 
NLO approximation leads to the estimate A^^*-* = 410 MeV^ that is indeed not far from the 
standard value 380 MeV (Appendix El). 

To determine 5^, we first fix the parameter A^ by employing the "conservative estimate" 
Ag(/i^ = 1 GeV^) = 0.4 GeV^. In QCD the value of this parameter was estimated in the QCD 
SR approach 22| and also using lattice data 



A 



(g(O)V^g(O)) ^"^^^"'"^ (g(0) {iga.^Cnqm 



mm) 



limit 



2(g(0)g(0)) 



0.4 ±0.1 GeV^ [20, 22 
0.4 - 0.5 GeV^ [6j 



(A.i; 



A brief review of the different estimates of A^ is given in 
The evaluation of the SR for the quantity r, 



1 - ^ ■ exp {-ml,/M'^) 



t2\ nf2 



27v^ (G^) 2aJl GeV' . 
9 A2M2 ^2 L 0/ 



(A.2) 



for the standard value of the gluon condensate {G ) = {asGG)/n = 0.012 GeV , j2^, and 




FIG. 4: The solid line denotes r(M^,p, s) for the best-fitted parameters g = 1.1, sq ~ 3.24. The short- 
dashed line corresponds to the value r = 1.07, whereas the long-dashed line corresponds to r calculated with 
A3 = 200 MeV. 



Using the values = 265 MeV and A4 ^° = 340 MeV given in '21'| for Nj = 4, we re-calculated the values 
for Nf = 3, i.e., = 312 MeV and A^^o ^ ^gv. 
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with the fitting parameters, 
So = 3.24 GeV^ yields 



i.e., the coupling to the it', q 1.1, and the duality interval 



1.07(±0.01) ^5'' ^ 0.19 GeV for = 0.312 MeV. 



(A.3) 

The stability of the SR ()A.2|) with respect to the Borel parameter is rather good, as 
the solid line in Fig. |3] clearly effects. 

Note that adopting the popular option A3*^^°'' = 200 MeV, it simply imitates the Q;f'"°(/i^) 
behavior at intermediate scales /x^, providing r = 1.02(±0.01) and 5^ ~ 0.20 GeV^. This 
old estimate ji^ is also presented here for comparison (cf. long-dashed line in Fig. 0}. The 
parameter r is also rather sensitive to the value of (G^). For example, the new estimate of 
the mean value (G^) = 0.009 GeV"^, suggested quite recently in 2J|, leads to r ~ 1.1 and 
5"^ ~ 0.18 GeV^. Taking into account the last estimate, we derive the uncertainty A^^^ ~ 
0.01 GeV^. The other uncertainties, inherent in the SR method ()A.2|) . result to an overall 
effect of the same order: viz., ~ 0.01 GeVl Finally, we can establish for 5 an accuracy: 
52 = (0.19 ± 0.02) GeV^ for X] = 0.4 GeV^ 

In the same way we obtain from ()A.1|) the corresponding central values and uncertainties 
of 5^ for the higher value of = 0.5 GeV^ (5^ = 0.235 ± 0.025) and analogously for the trial 
0.6 GeV^ (52 = 0.29 ± 0.03 GeV^), the latter being of interest due to instanton 



value 



models 



The one-loop anomalous dimension of 5^ is 7t4 
other hand, the one-loop scale dependence of 5^(/U^) 



= 32/9 (see, 
is given by 



for instance, [13 )• On the 



(A.4) 



APPENDIX B: CZ DA NORMALIZATION POINT 

The authors of 01 used as a normalization point for the CZ DA the scale fi^ = 0.5 GeV^, 
which is significantly larger than the original one used by Chernyak and Zhitnitsky: /i^^ = 
0.25 GeV^. This latter and rather low normalization point is due to the fact that CZ employed 
the description of charmonium decays, where the characteristic virtuality of the pion is indeed 
of this low order as yU^z- Furthermore, in order to construct their model at such a low scale, 
they evolved the 2nd moments, determined at a scale of 1.5 GeV^, down to this scale using 
1-loop evolution equations with Aqcd = 100 MeV. The result is the well-known CZ DA: 



1 + ^Cf (2x-l) 



30x(l-a;)(l-2x)^ (B.l) 



However, if one wants to know the shape of this model at another scale, one has to evolve it 
to that scale. But a natural question arises: what evolution equation should be used to do 
that? 

From our point of view, the best solution would be to determine once and for all the value of 
the second Gegenbauer coefficient of the CZ DA at the standard QCD SR scale, /ig = 1 GeV^, 
which is, also according to the CZ arguments, rather close to the scale of the second moment, 
1.5 GeV^. In this sort of determination, one needs to evolve from the CZ scale /i^g to the 
scale = 1.5 GeV^ of the second moment (^^)cz using the same 1-loop evolution equations 
(with Aqcd = 100 MeV) as in their original paper This produces 
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But after restoring this way the CZ model at a scale of 1.5 GeV^, we should use for the 
evolution to 1 GeV^ the actual value of the 1-loop QCD scale, i.e., 312 MeV. This gives 



a^^{l^^ = 1 GeV") = 0.56. 



(B.3) 



APPENDIX C: RADIATIVE CORRECTIONS 



1. Structure of the NLO amplitude Ti(Q^, g^; /u , x 



Here we present diagram per diagram the results of the NLO Ti calculation, performed in 
in the Feynman gauge using the "naive-75 scheme" jsi |^. To make the presentation 



TABLE IV: NLO results for individual diagrams^ 



Diagram 



Ti = iVr[ log-part + Ci 



In 



-part 



Ci{x,io) 



X P 



X P 



© CROSS 



In 



+ In 



Co{u,Uj) (g> Va{u,x) 
u 



Co{x,uj) 



X P 



b) 



X P 

M.C. © CROSS 



In 



2 In 



Co{u,uj) <S) Vh{u,x) 

u 



Cq{u,uj) — Co{u,u!) (2)Vb{u,x)^— 2Co{x,u!) 

J u 

Vb{u,x)+ + K{u,x) 



Co{x,u;) 



+ {ln[l+a;( 
-I- (x ^ x)| 



u — u 



u Q'^[l + uj{x - x)] 



X P 



X P 



In 



Co{x,u;) 



Cq{x,uj) - Co{x,Lu) 



M.C. near a diagram denotes the mirror conjugated diagram, cross - the diagram with crossed quark legs. 
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more compact, the average virtuality Q and the asymmetry parameter u 

have been used, employing also the notation x = 1 — x. 

The results are expressed in terms of the LO coefficient function Cq{x,uj) (see Eq. (j2.4j) ) 
and its logarithmic modification Cq{x,u!), naturally appearing in NLO calculations, 

^/ N 1 1 / N 1 \n[l+u(x-x)] . , ^, 

Co X, uj) = ^-—^ r + ix^x); x, u) = ^ \ \ + (x ^ x), C.2 

I + uj[x — X) 1 + uj[x — X) 

and their convolutions with Va and Vb, the latter being parts of the LO ERBL kernel, 

y y 0(x — y) 

Vaiy^x) = 2Cf6'(x - y)- + {x ^ x,y y) ; Vb{y,x) = 2Cf h {x x,y ^ y) ; 

X X X — y 



Vo{y, x) = Va{y, x)+ + Vb{y, x)- 



2Cf 



9{x-y)^[l + 



X 



x-y 



+ {x ^ x,y ^ y) 



. (C.3) 



Here, the +-form of a distribution V{y, x) is defined in the common way: 

V{y, x)-^- = V{y, x) — 5{y — x) / V{u,x)du. 



(C.4) 



The expression for diagram b) requires a more complicated construction involving K: 



1 

K{y,x) = CpOix-y)- 

X 



{x ^ x,y y) . 



(C.5) 



The log terms, containing the ultraviolet scale /i^ in Table IIVI are completely cancelled out 
on account of additional diagrams with self-energy corrections to the quark legs j^, with 
contributions of the form ln(yu|.//i|.) ■ Co(x,ci;). The cancellation of the ln(/i|^) terms for the 
full set of diagrams is a consequence of the Ward identity in QED. Finally, collecting the log 
terms from all diagrams, one obtains ln(Q^//i|') ■ Cq ® Vq in accordance with Eq. ()2.5|) . 

To perform the (formal) limit — ^ in ()2.5|) . one has to take a; — >^ 1 in the formulas of 
Table IIVI giving rise to the known expression for Ti |^ , [3| '■ 



T{Q^,0;fil;x) 



Nt 



X An 



In 



/ip J U u 



^0(^,0;) + 



X 



In X 



X Inx 



- 9 



X 



+ (x-.x) + 0(i,); 



Cl(g^O;x) 



1 2 a;lnx 

In X 9 

1 — X 



+ (x x). 



At the scale /i| = = Q^, this leads to 



(C.6) 
(C.7) 

(C.8) 

72 ' ' 900 ' ^'^''^^ 
where a„ = a^[Q'^^ and /i„ gives the size of the leading radiative corrections to the contri- 
bution of the nth Gegenbauer eigenfunctions '?/'„(x), entering the expansion of ^-nix) (see Eq. 



1 



1 02 + 04 + 



«s(Q')Cf 



47r 



-5 + 02/^2 + 04/^4) 



W)80- 

g2 27. 



K = ^Q^Ci{Q\{}]x)®i)n{x) 

O X 



295 



hi 



10487 
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2. NLO coupling constant 

Let us start with the RG equation for the rescaled running couphng as (Q^) 

das (Q^^ 



P{asiQ'),Nf) 



rfln(Q2) 

where Nj is the number of active flavors and the modified /3-function reads 



(C.IO) 



with Cfc = Cfc [Nf) 



P{as,Nf 

bo m 

bo (iV/) 



-al{l + cias + C2al + ...) , (C.ll) 
, and the standard /3-function coefficients are given by 



llNc-2Nf _ 



(C.12) 
(C.13) 
(C.14) 

Following here SY, we use the strong coupling constant in Sec. Illll in a "Particle Data Group" 
(PDG) form, which is the expanded second-order iteration of the 2- loop equation ()G.10|) : 



bi (Nf) = yiV,2 - (2Cp + Hat,^ n 



2857 5033 325 

Nf H N^f 

2 18 54 ^ 



a 



PDG 



with 



{Q^Nf) 



An 



In 



60 iNf)LiQ^ 



;Li(g2)=ciln [L{Q^)]-c^^ = cl-C2. 



K^NfY 

where we fixed A5 = A {Nj = 5) with the help of 

as {Q^ = (91.2 GeV)2, A^^ = 5) = 0.118 . 



(C.15) 



(C.16) 



(C.17) 



Matching this coupling at the Nf = 4 threshold, Q4 = 10 GeV, and analogously at the Nf = 3 
threshold, Q3 = 2.4 GeV, we arrive at 



A:, = 380 MeV . 



(C.18) 



But one can (we actually did this already in Sec. I1V|) use instead the exact solution of the 
2-loop RG equation, rather than the PDG-booklet expression. This exact solution can be 
expressed in terms of the quantity to read 



1 



a.(g2) 



Ci In 



1 



a,(Q2 



+ Ci 



HQ')- 



(C.19) 



As has been shown in |2^ the two-loop running coupling in QCD, being the solution of this 
equation, can be written via the Lambert W function: 

-A\Nf)-""^^^'' 



a 



2-loop 



An 



bo (Nf) ci 



1 + Wi 

Cl 



l/ci' 



(C.20) 
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The difference between tlie PDG form and tlie exact solution varies from 3.5% at = 10 GeV^ 
to 18% at = 1 GeV^ wlien one uses tlie same value of A = 380 MeV for both functions. In 
a real situation, the values of A in the two cases are different and are fixed in accordance with 
the standard boundary condition (irrTTfl . so that A™^ = 380 MeV and Ag"'""^ = 408 MeV. 
The deviation between the two forms becomes less pronounced at higher Q"^ and varies from 
0.7% at = 10 GeV^ to 10% at = i Q^y^ _ pig_ ^ 



0.45 



. 4 



0.35 



0.3 



0.25 




FIG. 5: Q^-dependence of two different 2-loop running couplings, af'^°°P{Q^) and a™^((5^): the solid 
line corresponds to the exact solution of the RG-equation, Eq. IIC.20|) : the dashed line to the PDG-booklet 
approximate form, Eq. (|C.15|I . 



It should be realized from this comparison that the PDG formula ()(I15j) is afflicted by a 
large error at = 1 GeV^ and for that reason it is preferable in the low Q^-region (< 2 GeV^) 
to use the exact formula ()G.20|) . It is worth to note here that in the case of treating the heavy- 
quark thresholds in a more accurate way, as done in [i^, the deviation between the PDG 
expression and the exact formula for the strong coupling becomes lower - at = 1 GeV^ 
about 5%. 



APPENDIX D: THE NLO EVOLUTION OF DA 

The ERBL evolution equation and its kernel have, respectively, the following structure (for 
more details, see, for example, |3^ ) 



In /i^ 

V{x,u;as) = [-^)^o{x,y) + (^^^ Vi{x,y) + ... 

Its eigenvalues 7„(a;s) and eigenfunctions \l'„(a;; a^) are defined through 

V{x,u;as) ^"^niu; as) = jnias)"^ n{x] as) ■ 



We use the following notations for its eigenvalues |31l 



ln{as) 



4:71 



fa \ *^ 

7oH+ \^) 7iH + --- 



(D.l) 
(D.2) 

(D.3) 
(D.4) 
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where the sign "— " just allows one to work with positive numbers 70 and 71 (n)'^: 



70 (0) 
7o(2) 

7o(4) 



0, 

100 

"9"' 

728 
I5"' 



7i(0) 
7i(2) 

71 (4) 



0, 

34450 830 



N 



243 81 
662846 31132 



3375 



2025 



Nf 



(D.5) 
(D.6) 

(D.7) 



The two-loop eigenfunctions ^ 2-i°°P(x;a,) of the ERBL equation (IIXT| can be expanded in 
terms of the one- loop eigenfunctions ipni^) (see Eq. (j2.8j) ). The approximate solution (j2.6j) in 
NLO is then given by 



^a„(/io)E(n,/i^ 



4tt 



k>n 



exp 



as ifJ-o) 



7„(47ra/6o) da 



(D.8) 
(D.9) 



with ao = 1, see, e.g., j32]. The "diagonal" part (in the {^n} basis) of ()D.8|) . expressed by 
the standard RG exponent, is the exact part of this solution, while the "non-diagonal" part 
is taken in the NLO approximation. The coefficients (i„fc(/i^, /Iq), corresponding to the non- 
diagonal part, fix the mixing of the higher harmonics, k > n, due to the fact that the matrix 
of the anomalous dimensions is triangular in the {ipn} basis. The exponent E{n,n'^) in ()D.8|1 
can be written explicitly (the scale /Iq can be fixed at some arbitrary value ~ 1 GeV^) to read 





7o{") 


'1 + cia,(Q2)" 






1 + cia^(/ig)_ 



uj(n) 



7i(n)&o -7o(ra)&i 



(D.IO) 
(D.ll) 



Evolving according to ()D.8|) . the Gegenbauer coefficients 02, change from the scale /Xq = /i^y^ 
to the scale as follows 



a2 ^ f/2(Q', a2) = a2^(2, Q^) + -J^d^^{Q\ iil) , 



04 U^iQ"^, 02, 04) = a4i?(4, Q"^) + 



47r 

asm 



+ a2E{2,Q')d2i{Q\fil 
Here the NLO mixing coefficients are (A; = 2, 4 > n = 0, 2 



[o^o4(Q^/io) 







r«s(Q')i 


[7o(fc)-7o(n)]/(26o) 


lo(,k) - 7o(n) - 2bo 




."s(Ai^). 





(D.12) 



(D.13) 



(D.14) 



^ The overall factor 1/2 is due to historical reasons in Ea. (|D.4|) . it is absent in the notations of the original 

paper [s^, whereas the factor (—1/2) is absent in MuUer notations [ssl . 
® Let us remind the reader that we use the values of 02 and 04 fixed at the scale /igy — 5.76 GeV^ as an input; 

this is done in order to facilitate comparison with the SY results. 
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where the values of the first few elements of the matrix Mnk are 

Mq2 = -11.2 + 1.73Nf, Mo4 = -1.41 + 0.565A^/, M24 = -22.0 + 1.65A^/ . (D.15) 

Analytic expressions for Mnk have been obtained in js^]- Using them one can estimate that 
the accuracy of ()D.15|) is of the order of 1%. 

One appreciates from Eq. ()D.8|) that the NLO evolution inevitably generates higher har- 
monics. Even in our case, where we have as a starting point only two harmonics 02 (/ig) 7^ 
and a4(/iQ) 7^ 0, the evolution to the scale produces a^i^Q^) 7^ for all k > 6. As one can 
see from Eq. ()D.8|) . for k > 6 these harmonics are of NLO {ak{Q'^) ~ as(Q^)). For this reason 
and owing to the enormous computional efforts needed for this task, we take into account only 
the complete NLO evolution of the first two nontrivial harmonics. 



APPENDIX E: EXPRESSIONS FOR THE NLO LCSR TRANSITION FORM FAC- 
TOR 

We employ a similar formalism as that used in 0, However, to make the present inves- 
tigation self-contained, we provide explicit expressions, given also that the formulas provided 
in are actually incomplete and only partially contained in All in all, the NLO LCSR 
transition form factor is 



F^^*^{M\Q\a2,ai) = F^^o^(M^ Q^ as, 04) + 



271 



(E.l) 



with 



G^'°(Q^ M') + t/2(g^ a2)Gf )^°(Q^ M') 



' ^Tw-4/,n2 Q\„—T7i-JC I ^ i_rTw-4/^2N 



m 



p JO 



p'-^-^Q^ S)e^dS + ^H'-^-^Q 

Q 



(E.2) 



where the evolution functions U2{Q^,a2) and ?74((5^, 02, ^4) are described in Appendix iDl We 
also define the following LCSR functions 



g2^0RDER^g2^^2 



n2 /-so rnf^^S 

^ / „ORDER/'/n2 C\ „ .,■> J C I UORDERY/n2N 



^0RDER^g2^ 



^ORDER^gi^ S)e^dS + if°^°™(Q^ 



so 
00 



^0RDER^,n2 n\Q_jl^ 

s 

Tw-4/,n2 ci\Q 



s 



(E.3) 
(E.4) 
(E.5) 
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for /c = 0, 2, 4 and set order= LO or NLO. The spectral densities in LO are 



6ms 



{U + S)3 ' 

36us(u^ — 3ms + s^) 

90ms(m^ - IOm^s + 20m2s2 - IOms^ + s^) 



(u + s 



|7 



and in NLO (see 0) 



—us 



u + s)'^ 
—s 



15 - TT^ + 3 In 



u 



4(m + s) 



10(m + s)7 



(E.6) 
(E.7) 
(E.8) 

(E.9) 

+ 12(s^-3sM + M^)Mln- (25 + 61n-) ; (E.IO) 
91s^ - 2(54 1 3 + 757r^)s^M + 125(541 + \2-n'^)s^u^ 

- 100(901 + 307r2)s\^ + 150(193 + 107r2)sM^ 

- 15(109 + lOvr^)^^ + 30(s^ - IOs^m + 20s V - IOsm^ + u^)u 

(E.ll) 



25s^ - 8(95 + 37r^)s^M + 36(25 + 27r^)sM^ - 12(5 + 2-n'^)u^ 



X In - f91 + 151n- 



Contributions from higher states in LO are given by 

37/2 



(u + So)2 ' 
3^2(^2 -8uso + 6s2) 

3u'^{u^ - 2Au^so + 90u'^sl - SOus^ + 15s^ 
{u + So)6 



(E.12) 
(E.13) 
(E.14) 



and in NLO: 



-15^2 



3so So _^ 3 (2m + So) So 



- 3 



2(m + so)2 (m + so) u 2(m + so)2 
So 



2 

M / 3 



In^V-" 



Li2 



u 



5u 



u 



+ 

+ 



48(m + so)4 
3so 



1-ln^) ln(l + ^^ 

u / u' 



(59m^ - 352^^50 + 564Msg - 72so^) 



(E.15) 



2(m + So)^ 
1 

3 



In^r-!^ 



(12^3 + 4mso' + 50=^) 

(5u^ + 4so (47m3 - 3so (18^^ + 3uso + So^))) In 

5 



So 

u 



Li,f-f5-)-(i-i„*)Mi + £i)l+2,„(i + ^); 

* "■ ' V U / M J 4 So 



u 



(E.16) 



In the last part of this exposition, we use u instead of in order to make the formulas more compact. 



23 



H, 



NLO, 



U 



U 



600(u + So) 



10487^^ - 149418m^so + 678285m^so^ - 888520m^so^ 



+ 259260mso^ - 6300so^ 



+ 



3s 







2(n + So)6 
u 

+ — — : (44m^ + 2124mSo - 7575^^50^ + 8000^^50 



In 



vr 



10(m + So) 



- 1425uso^ + 30so^) In — - 3 fin 
^ u \ 



-^(-7)-(-'"?)'"(i- 



u 

So 

u 



u 

5 So 



The twist-four spectral density is 



P 



u,s) 



806^{u) 2u^s{s-u) 



3u {u + sY ' 
where 5'^{u) is taken from Eq. ()A.4|) with /x^ = u. This produces 

Tw-4/ N _ 80^^(m) m^(2so -m) 



(E.17) 



(E.18) 



(E.19) 
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